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Bats are connected with the increasing numbers of emerging and re-emerging viruses that may break the species barrier and 
spread into the human population. Coronaviruses are one of the most common viruses discovered in bats, which were considered 
as the natural source of recent human-susceptible coronaviruses, 1.e. SARS-COV and MERS-CoV. Our previous study reported 
the discovery of a bat-derived putative cross-family recombinant coronavirus with a reovirus gene p/0, named as Ro-BatCoV 
GCCDC1. In this report, through a two-year follow-up of a special bat population in one specific cave of south China, we 
illustrate that Ro-BatCoV GCCDC1 persistently circulates among bats. Notably, through the longitudinal observation, we 
identified the dynamic evolution of Ro-BatCoV GCCDC1 in bats represented by continuously recombination events. Our study 
provides the first glimpse of the virus evolution in one longitudinally observed bat population cohort and underlines the 
surveillance and pre-warning of potential interspecies transmittable viruses in bats. 
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Bats (Order Chiroptera) are thepg only group of mammals 
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e., nocturnal animals (Calisher et al., 2006). Bats have been 
the focus of particular interest not just because of their 
ecological interest, but also because of their medical im- 
portance in public health, as they play significant role in 
evolving and transmission of zoonotic viral agents that 
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stance a substantial danger to human and animal well-beings. 
Bats possess special peculiarities and features which make 
them host to different viruses, such features including 
structural diversity and complex lifestyle, diverse mode of 
feeding ranging from insectivores (feed on insects), frugi- 
vore (fruit eater), nectarivore (sucking nectar), palynivore 
(eating pollen), piscivore (eating fish) while others feed on 
arthropods, vertebrates or practice haematophagy (feed on 
blood). Bats are social animals that live in a densely popu- 
lated roost, demographic and spatial structure of population 
with long life, which may significantly contribute to main- 
taining viruses as well as the high possibility of intra- and 
inter-species transmission of viral infections (Calisher et al., 
2006). In addition, excreta during migration may be an im- 
portant factor in the spread of viruses (Lopez-Roig et al., 
2014). 

Bats are connected with the increasing numbers of emer- 
ging and re-emerging viruses including Hendra virus (HeV), 
Nipah virus (NiV) and Ebola virus (EBOV), and also cor- 
onaviruses, etc. (Li et al., 2016). The studies of bats as re- 
servoirs for a number of mammalian coronaviruses came to 
the spotlight since the discovery of severe acute respiratory 
syndrome (SARS)-like coronaviruses in 2005 isolated from 
south China (Lau et al., 2005; Li et al., 2005), and subse- 
quently, a SARS-CoV receptor ACE2-using coronavirus was 
isolated from bat (Ge et al., 2013). Sister to SARS-CoV was 
the middle east respiratory syndrome-associated coronavirus 
(MERS-CoV) emerged in 2012, the close-related viruses of 
which was also found in bats (Lelli et al., 2013; Zumla et al., 
2016). Actually, a number of bat coronaviruses have been 
documented in recent years (Hu et al., 2017), and the MERS- 
CoV phylogenetically closely related bat coronaviruses 
(BatCoVs) HKU4 can also recognize human CD26 for cell 
entry, indicating origin of human-susceptible coronaviruses 
from bats (Wang et al., 2014). 

Coronavirus encrypts a non-segmented, largest positive 
sense single-stranded RNA genome known ranging from 27 
to 32 kb in length (Lai, 1996; Sokolova et al., 1996; Du et al., 
2016). In the subfamily of Coronaviridae, Alphacoronavirus 
and Betacoronavirus infect and cause diseases in mammals, 
while Gammacoronavirus and Deltacoroanvirus are mostly 
confined to avian infections and diseases (Jackwood et al., 
2012; Woo et al., 2012). Report has shown that coronaviruses 
possess low proofreading capability that can also exhibit 
comparatively high mutation and recombinant rates, which 
make them one of the most diverse, genetically distinct, and 
recently emerging groups of viruses (Minskaia et al., 2006). 
The continuous mutation and recombination of the cor- 
onaviruses probably bridge the species barrier and enable 
them to spread into the human population (Baric et al., 1995). 

Recently, we discovered a novel Rousettus bat coronavirus 
GCCDC1 (Ro-BatCoV GCCDC1) with full genome se- 
quenced, which exhibits a putative heterologous re- 
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combination event between a bat coronavirus and a bat 
orthoreovirus (Huang et al., 2016). We found that the p/0 
gene of bat orthoreovirus in Ro-BatCoV GCCDC1 genome 
has no homologous gene from any parental strain of cor- 
onaviruses. However, there is dearth of information whether 
this virus with the recombinant gene can persist and evolve 
in their hosts. 

Herein we carried out a longitudinal survey on one po- 
pulation of Roussetus leschenaulti bats in Yunnan Province, 
Mainland China, between 2014 and 2015. The persistence of 
the Ro-BatCoV GCCDCI1 in the bat population and its po- 
tential evolution was investigated. Our study provides the 
first glimpse of the virus evolution in one prospectively and 
longitudinally observed bat population. 


RESULTS 


The persistent prevalence of the Ro-BatCoV GCCDCI1 
in the bats 


A total of 568 rectal swab samples from R. leschenaulti bats 
were collected between May 2014 (n=118), October 2014 
(n=270) and October 2015 (n=180) in one cave located in 
Xishuangbanna, Yunnan Province, China (Table 1). The anal 
swab samples were screened individually for the presence of 
coronavirus using RT-PCR detection method. Coronavirus 
positive ratios were 39.8% (47/118), 38.8% (70/270) and 
35.6% (64/180) in bats sampled in May 2014, October 2014 
and October 2015, respectively (Table 1). These results 
showed a persistent prevalence of Ro-BatCoV GCCDC1 in 
this bat population. 


The evolution of RdRp from Ro-Bat GCCDC1 


Based on the identified PCR CoV-positive bat samples, the 
PCR amplicons with the length of 228 bp of Ro-BatCoV 
GCCDC1 RdRp were sequenced, which showed the con- 
servation of this newly identified bat CoV. We analyzed the 
RdRp sequences which had the corresponding p/0 sequences 
available, thus we have 97 pairs of RdRp and p/0 from 97 bat 
specimens. The phylogenetic analysis of the RdRp sequences 
revealed two groups of Ro-BatCoV GCCDC1 designated 
Group I (RdRp-I) and Group I (RdRp-II) (Figure 1A). The 
result from our study revealed that all the Ro-BatCoV 
GCCDC1 RdRp from the first round sampling in May 2014 
belonged to the RdRp-I, while subsequent sampling in Oc- 
tober 2014 and 2015, respectively, showed the involvement 
of viruses featured with RdRp-IlI. 


The evolution of p10 gene in Ro-BatCoV GCCDCI1 


Based on the result of PCR amplification using specific 
primers, the nucleic acid sequences of p/0 genes (420 bp) of 
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Table 1 Summary of sampling detail and coronavirus infections in bats” 
Bats with swab 


Bats with CoV 
samples collected 


Sampling time positive (%) 


May 2014 118 47 (39.8%) 
October 2014 270 70 (38.8%) 
October 2015 180 64 (35.6%) 


a) *, All the positive samples were defined based on the PCR amplicons 
with the length of 228 bp of Ro-BatCoV GCCDC1 RdRp as described in 
Table 2. 97 RdRp sequences out of the 181 samples which had the corre- 
sponding p/0 sequences were utilized for the GCCDC evolution. 
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Figure 1 Phylogenetic analyses of Ro-BatCoV GCCDC1 for three 
sampling periods. Evolutionary trees of RdRp (A) and P/0 (B) constructed 
for the samples collected from May, 2014 to Oct., 2015 has indicated by 
different colors as labeled. The maximum-likelihood trees were constructed 
with 1,000 replication available in MGEA6. 


Ro-BatCoV GCCDC1 were subjected to phylogenetic ana- 
lysis, which displayed high level of sequence similarity to 
Ro-BatCoV GCCDC1 p/0 genome segments (Huang et al., 
2016). In pairwise comparison of nucleotide sequence, the 
nucleotide sequence identity was 99% compared to Ro- 
BatCoV GCCDC1. However, based on the phylogenetic 
analysis of all the available Ro-BatCoV GCCDC1 in public 
domain, the p/0 genes were divided into two groups: Group 
(A) and Group (B) (Figure 1B). All the three times of sam- 
pling Ro-BatCoV GCCDC1 covered both Group A and 
Group B, which indicated the persistence of the groups of 
p10 in Ro-BatCoV GCCDC1. 

When compared to p/0 genes from traditional orthor- 
eoviruses (Figure 2), the p/0 from the typical six Ro-BatCoV 
GCCDC1: GCCDC1/2-53 and GCCDC1/2-110 sampled in 
October 2014, and GCCDC1/3-19 and GCCDC1/3-31 in 
October 2015 and our previously identified Ro-BatCoV 
GCCDC1 strains 346-KU762337 and Ro-BatCoV GCCDC1 
strains 356-KU762338, are clustered into one independent 
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branch indicating the uniqueness of the bat coronavirus-de- 
rived p10 gene. However, we can still observe that p/0 genes 
in 2-53 and 3-31 belong to Group A and p/0 of other four 
samples are in Group B. 


The potential dynamic recombination of Ro-BatCoV 
GCCDC1 


To further investigate the dynamic evolution of Ro-BatCoV 
GCCDC1 in the longitudinally collected bat samples, we 
further analyzed the potential recombination of the viruses 
by constructing a model that pairs of RdRp and p/0 genes, 
each pair of which is sequenced from single Ro-BatCoV 
GCCDC1 sample. Five different nonsynonymous mutations 
were found in all available RdRp protein segments (Figure 
3A). However, only two mutations were featured by the two 
groups of RdRp (T-648 and I-666 for RdRp-I and A-648 and 
M-666 for RdRp-IJ). Comparison of the translated amino 
acid sequences of p/0 gene of both Group A and B of Ro- 
BatCoV GCCDC1 identified four nonsynonymous muta- 
tions at position 7, 14, 19 and 34, respectively (Figure 3B). 
Group A p10 was featured with the combination of M-7, T- 
14, A-19, and I-34, while Group B possesses substitutions V- 
7, T-14, N-19 and S-34 in the four positions. When we paired 
the corresponding RdRp and p/0 genes groups with their 
featured substitutions, there were four different combina- 
tions of the putative assembled genomes in all the bat sam- 
ples tested (Figure 4A). Ro-BatCoV GCCDCl1 virus IA 
possesses Group I of RdRp gene and Group A of p/0 gene. 
Ro-BatCoV GCCDC1 virus IB is featured with Group I of 
RdRp gene and Group B of p/0 gene. Ro-BatCoV GCCDC1 
virus IIA has Group II of RdRp gene and Group A of p/0 
gene. Ro-BatCoV GCCDC1 virus IIB characterizes by 
Group II of RdRp gene and Group B of p/0 gene. During the 
first round of sampling in May 2014, only viruses IA and IB 
were discovered in the bat population (Figure 4B). However, 
five months later in October 2014 IIA and IIB recombinants 
can be found in the bat swabs. Further, the third round 
sampling in October 2015 showed that while the four com- 
binations were persistent in the bat population, the two 
viruses ITA and IIB became more dominant. 


DISCUSSION 


In our previous study, we demonstrated that a novel cor- 
onavirus Ro-BatCoV GCCDC1 encoded a putatively re- 
combinant nonstructural protein p10 derived from reovirus, 
which induces the formation of syncytia in infected cells 
(Huang et al., 2016). Sequel to this, we now report that Ro- 
BatCoV GCCDC1 is continuously circulating in bats from 
Yunnan Province in China. We observed high prevalence of 
the Ro-BatCoV GCCDC1 on a consistent trend for a period 
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Figure 2 Phylogenetic relationships of Ro-BatCoV GCCDC1 compared to other viruses. Phylogenetic trees of the p/0 sequence of Ro-BatCoV GCCDC1 
together with the p/0 genes from other reoviruses. The phylograms were constructed by using MEGA6 software package to determine the maximum- 


likelihood tree for the sequence. 
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Figure 3 Multiple sequence alignment of Ro-BatCoV GCCDC1 genes based on 


eS Se OO Oe Oe Oe 


ui SVLCUOG!TVKPUI SSDS 


4 ald phd ald 2) 4 ed old 
ALATA CATA CA CATE CAG 


ee ee 


SSe-SeSeeeeree — 
AAS NS>NSS>Soe> 


Il 


eh he rhe she ole he old hd od phe odd 2h dp 


( G « 
( (r ( 
( (y « 
( G ¢ 
CGGI 
.CGOGGEI 
CGGG 
CGGG 
( G ¢ 
( G ¢ 
( Gi 
( (; 
( G 


we ee eee ee 


DOLD oOHaDoHL LA 


NIAC WNW 


A TROGIR CIEE Y CC YIEK F 

TVGLLLLY CLEYCCYLAFKASO! Ctr rion PRET rj 
IVGLLLLYVCLFEYCCYLKFKASQAKKTYRKELIS RQLYAPPREISHY 
IVGI VCLFYCCYLKFKASQAKKTYRKELIS - RQIYAP RETSHY 
IVGI VCLFYCCYLKFKASQAKKTYRKELISLTTRQIYAPPREISH\ 
IVGI VCLFYCCYLKFKASQAKKTYRKELISLTTRQIYAPPREISH\ 
IVGLLLEVCLEYCCYLKEKASQAKKTYRKELISLTTRQIYAPPRETSHY | (ay 
IVGLLLLVCLFYCCYLKFKASQAKKTYRKELISLTTRQIYAPPREISH\ 
IVGLLLLVCLFYCCYLKFKASQAKKTYRKELISLTTRQIYAPPREISHV 
IVGI VCLFYCCYLKFKASQAKKTYRKELISLTTRQIYAPPREISHY 
IVGI LVCLFYCCYLKFKASQAKKTYRKELISLTTRQIYAPPREISH\ 
IVGI VCLFYCCYLKFKASQAKKTYRKELISLTTRQIYAPPREISHY 
IVGI VCLEYCCYLKEKASQAKKTYRKELISLTIRQLYAPPREISHY | (B) 
IVGI VCLFYCCYLKFKASQAKKTYRKELISLTTRQIYAPPREISHY 


samples from three sampling periods. Sequence alignment of deduced 


amino acid of RdRp (A) and p10 proteins (B). Consensus and sequence logo for conserved amino acid residues are shown on top of the alignment. The 


residues with color indicate the position of mutation. 


of 2 years that is relatively high in proportion to the popu- 
lations under study. The cross-family recombinant p/0 gene 
stably exists in the novel virus. Notably, we observed the 
dynamic evolution of this virus in the special bat population. 
Our study provides beneficial recommendation for the sur- 
veillance and pre-warning of potential interspecies trans- 
mittable viruses in bats. 


The recombination events in coronavirus involving inter- 
change of intra-family genetic materials has been reported 
(Baric et al., 1995). In our previous study, a traditional re- 
ovirus gene p/0 located at the upstream of N protein and 
downstream of NS7a of Ro-BatCoV GCCDC1 genome 
showing the ability of coronavirus to undergo a heterologous 
inter-family recombination (Huang et al., 2016). In our 
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Figure 4 Recombination analysis of Ro-BatCoV GCCDC1. A, The re- 
combination model for Ro-BatCoV GCCDC1 based on the different groups 
of RdRp (Group I and Group I) and p10 proteins (Group A and Group B). 
The groups are based on the different combinations of polymorphic amino 
acid residues in Ro-BatCoV GCCDC1. B, The ratios of different re- 
combinants of Ro-BatCoV GCCDC1 revealed two new groups, ITA and IIB 


current study, we observed a persistent prevalence Ro-Bat- 
CoV GCCDC1 in the bat population, and we also found that 
the recombinant p/0 is a constitutive gene in this virus. We 
did not get all the p/0 sequences from the RdRp positive 
samples. This may be due to a lower amplification efficiency 
of the p/0 primers. Thus, the deep sequencing of the whole 
genomes and also the virus isolation from the original sam- 
ples are still needed to elucidate more features of Ro-BatCoV 
GCCDC1. 

The emergence of recombinant p/0 gene in Ro-BatCoV 
GCCDC1 could be due to unhindered access of reoviruses in 
bat, being a prehistoric mammal which might have favoured 
viral integration and persistence. p/0 persistence may have 
utilized intrinsic factors mimicking ability to incorporate its 
genome into coronavirus or of the host’s cells. p/0 has been 
reported to share similarities with viroporins, which has the 
ability to form hydrophilic channels that allow low mole- 
cular weight hydrophilic molecules to cross the membrane 
thereby interrupt membrane integrity, breakdown ionic gra- 
dients and eventual release of essential compounds from the 
cell (Carrasco et al., 1989; Carrasco et al., 1993). This could 
lead to series of events that favor viral replication, encourage 
viral transmission and persistence of Ro-BatCoV GCCDC1 
in bat population. Interestingly, the phylogenetic analysis 
revealed that two distinct groups of p10 with featured amino 
acids persistently exist in Ro-BatCov GCCDC1. However, 
the comparison of the two p10 in the pathogenesis of bat 
corornavirus need further investigations. 

Summarily, the new Ro-BatCoV GCCDC1 virus is not 
only persistently prevalent in the bat population, but also 
exhibits continuous evolution. Probability also exists that 
Ro-BatCoV GCCDC1 might be under diverse selective 
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pressures. The evolving bat coronaviruses could have seized 
the advantage of conserved mammalian receptors for cellular 
replication and other biological pathways that go through 
radiation over a long time which could boost the capability 
for conveyance of bat-related viruses to other animals 
(Calisher et al., 2006). Our study benefits the understanding 
of dynamic virus evolution in special bat population and 
underlines the surveillance of potential interspecies viruses 
in bats. 


MATERIALS AND METHODS 


Sample collection, storage and transportation 


Anal swab samples were collected from Rousettus lesche- 
naulti (R. leschenaulti) bats captured on three occasions, 
May 2014, October 2014 and October 2015, in one cave 
located in Xishuangbanna, Yunnan Province, China. Rectal 
swabs were collected from each bat and placed in viral 
transport medium (VTM) in screw cap tubes, kept on cold ice 
(at 4°C to 10°C) and transported to the laboratory at local 
CDC where they were frozen in liquid nitrogen. The VTM 
containing Earle’s Balanced Salt Solution, 5% bovine albu- 
min, 50,000 ng mL vancomycin, 50,000 ng mL | amika- 
cin, 10,000 units mL ' nystatin as described in our previous 
study (Huang et al., 2016). The samples were later trans- 
ported to the laboratory in dry ice and stored at —80°C until 
further analysis. 


RNA extraction 


VTM containing swab samples were manipulated by cen- 
trifugation at 8,000 r min ' for 5 min at 4°C using a Thermal 
centrifuge. Total RNA was extracted from 140 uL of the 
supernatants using the RNeasy mini kit (Qiagen, Germany) 
according to the manufacturer’s protocol. The RNA was 
eluted in 60 uL AVE buffer and stored at —80°C for sub- 
sequent use. 


Viral detection by nested-PCR 


The extracted total RNA was transcribed to cDNA using 
reverse transcriptase enzyme (Thermo, USA) and then 
screened for the presence of coronavirus RNA using pan- 
coronavirus RT-PCR degenerate primers. The primers were 
designed from a highly conserved RdRp region of cor- 
onaviruses (Table 2). Theoretically, this pan-coronavirus RT- 
PCR primers should amplify a 299 bp fragment of the 
polymerase gene of all coronaviruses in the first round and 
228 bp in a semi-nested PCR for the second round. If the 
samples were RT-PCR positive to coronavirus screening, 
specific p/0 primers (Table 2) that covered the whole p/0 
gene were used to screen for the presence of reovirus p/0 
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Table 2 


The primers for coronaviridae gene identification and amplification 
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Primer name Target Sequence (5’—3’) G+C (%) Ty CC) Pete or me -ampiueauon 
gene fragments 
RdRp-OF RdRp TGTTATTGGAACCACGAAGTTY TAY GGNGGNTG 42.1 51.3 OF+OR=299 bp 
RdRp-lF RdRp GTTTTGTTCTTATGGGTTGGGATTAY CCNAARTGYGA 42.1 52.3 IF+OR=2728 bp 
RdRp-OR RdRp TAGTAGCATCTCCGCTGCTAGTNCCNCCNGGYTT 52.6 53.4 
CoV-p10-OF plo GAAAGGGAGACTAAGAAGA 42.1 50.8 
OF+OR=715 bp 
CoV-p10-OR plo ACAGGCTATAACTGTAAATG 35.0 49.2 
CoV-p10-IF plo GCAGATGAGAATGACAAACCCAA 43.5 56.0 
IF+IR=447 bp 
CoV-p10-IR plo GCAGCAAGCAAGCAGAGCAA 55.0 57.4 


gene in Ro-BatCoV GCCDC1. The amplified PCR products 
were purified following electrophoresis in 1% agarose gels 
that revealed the expected amplicon in agarose gels. 


Virus isolation 


RT-PCR positive samples for coronavirus were cultured in 
human epithelial colorectal adenocarcinoma cells (CaCo-2), 
human lung carcinoma cells (A549), human epithelial type 2 
HeLa derivative cells (HEp-2), baby hamster kidney (BHK- 
21) and rhesus monkey kidney cells (LLC-mk2). Cells were 
grown in Dulbecco’s modified Eagle’s minimum essential 
medium (DMEM) supplemented with 10% or 20% fetal 
bovine serum (FBS) and incubated at 37°C with 5% CO. 
The cell lines were inoculated with the supernatant of PCR 
positive samples in 24-well plates, incubated for 2 h for virus 
adsorption, 300 uL fresh cell culture medium with 2% FBS 
was added and incubated at 37°C in 5% CO, incubator. The 
plates were observed daily for viral cytopathic effect (CPE). 


Molecular and sequencing analysis 


The used reference sequences were downloaded from Gen- 
Bank. Multiple sequence alignment of nucleic acid se- 
quences was performed using ClustalW and the Geneious 
program (http://www.geneious.com). For RdRp and pl0 
phylogenetic analysis, maximum-likelihood trees were con- 
structed with MEGA 6 (http://www.megasoftware.net), and 
1,000 bootstrap replicates were run. 


Compliance and ethics 
of interest. 


The author(s) declare that they have no conflict 


Acknowledgements We appreciate the great efforts of Drs. Yongming 
Zhou, Honghua Wen, Huaxing Liu, who participated in the collection of the 
bat samples. This work was supported by the National Key Research and 
Development Program of China (2017YFC1200202), the National Natural 
Science Foundation of China (81290342, 1461168030), the Major Special 
Projects for Infectious Disease Research of China (2016ZX10004222-003), 
and China National Grand S&T Special Project (2014ZX10004-001-006). 


Joseph O. Obameso was supported by CAS-TWAS President's Fellowship of 


the University of Chinese Academy of Sciences (UCAS) and The World 
Academy of Sciences (TWAS). George F: Gao is a leading principle in- 
vestigator of the NSFC Innovative Research Group (81621091). The funders 
had no role in study design, data collection and analysis, decision to pub- 
lish, or preparation of the manuscript. 


Baric, R.S., Fu, K., Chen, W., and Yount, B. (1995). High recombination 
and mutation rates in mouse hepatitis virus suggest that coronaviruses 
may be potentially important emerging viruses. Adv Exp Med Biol 380, 
571-576. 

Calisher, C.H., Childs, J.E., Field, H.E., Holmes, K.V., and Schountz, T. 
(2006). Bats: important reservoir hosts of emerging viruses. Clin Mic- 
robiol Rev 19, 531-545. 

Carrasco, L., Otero, M.J., and Castrillo, J.L. (1989). Modification of me- 
mbrane permeability by animal viruses. Pharmacol Ther 40, 171-212. 

Carrasco, L., Perez, L., Irurzun, A., Lama, J., Martinez-Abarca, F., Rodr- 
iguez, P., Guinea, R., Castrillo, J.L., Sanz, M.A., and Ayala, M.J. (19- 
93). Regulation of Gene Expression in Animal Viruses. (New York: 
Plenum Publishing Corp.), pp. 283-305. 

Du, J., Yang, L., Ren, X., Zhang, J., Dong, J., Sun, L., Zhu, Y., Yang, F., 
Zhang, S., Wu, Z., and Jin, Q. (2016). Genetic diversity of corona- 
viruses in Miniopterus fuliginosus bats. Sci China Life Sci 59, 604-614. 

Ge, X.Y., Li, J.L., Yang, X.L., Chmura, A.A., Zhu, G., Epstein, J.H., Mazet, 
J.K., Hu, B., Zhang, W., Peng, C., Zhang, Y.J., Luo, C.M., Tan, B., 
Wang, N., Zhu, Y., Crameri, G., Zhang, S.Y., Wang, L.F., Daszak, P., 
and Shi, Z.L. (2013). Isolation and characterization of a bat SARS-like 
coronavirus that uses the ACE2 receptor. Nature 503, 535-538. 

Hu, B., Zeng, L.P., Yang, X.L., Ge, X.Y., Zhang, W., Li, B., Xie, J.Z., Shen, 
X.R., Zhang, Y.Z., Wang, N., Luo, D.S., Zheng, X.S., Wang, M.N., 
Daszak, P., Wang, L.F., Cui, J., and Shi, Z.L. (2017). Discovery of a rich 
gene pool of bat SARS-related coronaviruses provides new insights into 
the origin of SARS coronavirus. PLoS Pathog 13, e1006698. 

Huang, C., Liu, W.J., Xu, W., Jin, T., Zhao, Y., Song, J., Shi, Y., Ji, W., Jia, 
H., Zhou, Y., Wen, H., Zhao, H., Liu, H., Li, H., Wang, Q., Wu, Y., 
Wang, L., Liu, D., Liu, G., Yu, H., Holmes, E.C., Lu, L., and Gao, G.F. 
(2016). A bat-derived putative cross-family recombinant coronavirus 
with a reovirus gene. PLoS Pathog 12, e1005883. 

Jackwood, M.W., Hall, D., and Handel, A. (2012). Molecular evolution and 
emergence of avian gammacoronaviruses. Infect Genet Evol 12, 1305- 
1311, 

Lai, M.M.C. (1996). Recombination in large RNA viruses: coronaviruses. 
Semin Virol 7, 381-388. 

Lau, S.K.P., Woo, P.C.Y., Li, K.S.M., Huang, Y., Tsoi, H.W., Wong, B.H.L., 
Wong, S.S.Y., Leung, S.Y., Chan, K.H., and Yuen, K.Y. (2005). Severe 
acute respiratory syndrome coronavirus-like virus in Chinese horseshoe 
bats. Proc Natl Acad Sci USA 102, 14040-14045. 

Lelli, D., Papetti, A., Sabelli, C., Rosti, E., Moreno, A., and Boniotti, M.B. 
(2013). Detection of coronaviruses in bats of various species in italy. 
Viruses 5, 2679-2689. 

Li, W., Shi, Z., Yu, M., Ren, W., Smith, C., Epstein, J.H., Wang, H., 


Obameso, J.O., et al. 


Crameri, G., Hu, Z., Zhang, H., Zhang, J., McEachern, J., Field, H.., 
Daszak, P., Eaton, B.T., Zhang, S., and Wang, L.F. (2005). Bats are 
natural reservoirs of SARS-like coronaviruses. Science 310, 676-679. 

Li, Z., Liu, D., Ran, X., Liu, C., Guo, D., Hu, X., Tian, J., Zhang, X., Shao, 
Y., Liu, S., and Qu, L. (2016). Characterization and pathogenicity of a 
novel mammalian orthoreovirus from wild short-nosed fruit bats. Infect 
Genet Evol 43, 347-353. 

Lopez-Roig, M., Bourhy, H., Lavenir, R., and Serra-Cobo, J. (2014). Ser- 
oprevalence dynamics of european bat lyssavirus type 1 in a multi- 
species bat colony. Viruses 6, 3386-3399. 

Minskaia, E., Hertzig, T., Gorbalenya, A.E., Campanacci, V., Cambillau, 
C., Canard, B., and Ziebuhr, J. (2006). Discovery of an RNA virus 
3’—5' exoribonuclease that is critically involved in coronavirus RNA 
synthesis. Proc Natl Acad Sci USA 103, 5108-5113. 

Sokolova, T.M., Uryvaev, L.V., Selivanova, T.K., Bobrova, O.V., Lebedev, 
A., and Bystrov, N.S. (1996). Recombination of alpha-viruses during 
mixed infection in cultured cells. Formation of hybrid forms of RNA of 


Sci China Life Sci T 


Venezuelan equine encephalomyelitis virus and Karelian fever virus in 
the region of genes for nucleocapsid and envelope proteins. Vopr Vir- 
usol 41, 245-252. 

Wang, Q., Qi, J., Yuan, Y., Xuan, Y., Han, P., Wan, Y., Ji, W., Li, Y., Wu, Y., 
Wang, J., Iwamoto, A., Woo, P.C., Yuen, K.Y., Yan, J., Lu, G., and Gao, 
G.F. (2014). Bat origins of MERS-CoV supported by bat coronavirus 
HKU4 usage of human receptor CD26. Cell Host Microbe 16, 328-337. 

Woo, P.C.Y., Lau, S.K.P., Lam, C.S.F., Lau, C.C.Y., Tsang, A.K.L., Lau, J. 
H.N., Bai, R., Teng, J.L.L., Tsang, C.C.C., Wang, M., Zheng, B.J., 
Chan, K.H., and Yuen, K.Y. (2012). Discovery of seven novel mam- 
malian and avian coronaviruses in the genus deltacoronavirus supports 
bat coronaviruses as the gene source of alphacoronavirus and betacor- 
onavirus and avian coronaviruses as the gene source of gammacor- 
onavirus and deltacoronavirus. J Virol 86, 3995-4008. 

Zumla, A., Chan, J.F.W., Azhar, E.I., Hui, D.S.C., and Yuen, K.Y. (2016). 
Coronaviruses—drug discovery and therapeutic options. Nat Rev Drug 
Discov 15, 327-347. 


